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ABSTRACT 

A 3D algebraic grid generation scheme is presented for generating the grid points inside gas 
turbine combustors with inclined slots. The scheme is based on the 2D transfinite interpolation 
method. Since the scheme is a 2D approach, it is very efficient and can easily be extended to gas 
turbine combustors with either dilution hole or slot configurations. 

To demonstrate the feasibility and the usefulness of the technique, a numerical study of the 
quick-quench/lean-combustion (QQ/LC) zones of a staged turbine combustor is given. Preliminary 
results illustrate some of the major features of the flow and temperature fields in the QQ/LC zones. 
Formation of co- and counter-rotating bulk flow and sandwiched-ring- shape temperature fields can 
be observed clearly, and the resulting patterns are consistent with experimental observations typical 
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of the confined slanted jet- in-cross flow. Numerical solutions show the method to be an efficient 
and reliable tool for generating computational grids for analyzing gas turbine combustors with 
slanted slots. 
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INTRODUCTION 


The current issue of environmental pollution imposes an urgent need in the reduction of NOx 
emissions from gas turbine engines. One way to reduce such NOx emissions is to use the staged 
turbine combustor (STC) concept [1*4]. One STC combustor is schematically shown in Figure 1. 
It consists of a fuel nozzle (FN), a rich-bum (RB) zone, a converging pipe, a quick-quench (QQ) 
zone, a diverging pipe, and a lean-combustion (LC) zone. 

There are basically two approaches for modeling chemically reactive fluid flow and heat 
transfer phenomena in the STC combustor shown in Figure 1. The first is to consider the whole 
combustor as a single unit. The second is to isolate each zone and then connect them through 
inlet/outlet boundary conditions. This allows the zones to be treated differently. For example, the 
FN/RB zone can be assumed to be axisymmetric (with swirl), while the QQ/LC zones must be 
modeled in three-dimensional (3D) due to the jet-in-cross flow. From computational efficiency 
view point, the second approach is preferred and was used in our approach [5 -7] in which the FN/RB 
zone was calculated as an axisymmetric problem with a finer mesh to resolve the airblast nozzle 
passage flow fields [5,6] and the QQ/LC zones was modeled as a 3D problem using a coarser sector 
mesh to isolate one slot [5,7]. 

The purpose of this paper is to present an algebraic 3D grid generation scheme for gas turbine 
combustors with inclined slots. The method used to generate the grids is based on the 2D transfinite 
interpolation (TI) method [8-1 1]. Since the scheme is a 2D approach, it is very efficient and can 
easily be extended to gas turbine combustors with either dilution hole or slot configurations. 
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To demonstrate the feasibility and the usefulness of the approach just outlined, numerical 
solutions for the chemically reactive flow and heat transfer fields inside the QQ/LC zones of the 
STC are obtained using a modified version of the KTVA-II code [12]. Input conditions are chosen 
to be similar to those encountered in advanced combustion systems. 

DESCRIPTION OF PROBLEM 

We will focus on the QQ/LC region of the STC combustor shown in Figure 1. In this region, 
cool dilution air is injected into the QQ zone through slanted slots. In the model, to minimize the 
effect of the inter-zone boundary conditions, the inlet boundary is extended 3 inches upstream of 
the convergence section (see Figure 1 for the FN/RB zone and QQ/LC zones computational sections). 
Likewise, to avoid flow contamination due to outflow boundary conditions, the LC zone is extended 
6.6 inches downstream of the divergence section. Figure 2 shows the geometry of the QQ/LC zones. 
In this figure, D and L qujck are the diameter and the length of the QQ zone, respectively. 

There are equally- spaced inclined slots located around the perimeter of the QQ zone. The 
center of the slots is located at D/3 downstream of the inlet of the QQ zone. The slot opening area 
is determined such that the jet-to-mainstream mass flow rate ratio is a chosen constant with a 
predetermined slot orifice discharge coefficient. Due to geometric symmetry, only one slot is 
modeled resulting in a sector as shown in section AA of Figure 2, in which the sector angle is 45° 
for a 8 slots configuration. The slot is symmetrically located at the center plane of the sector. 

THE GRID GENERATION SCHEME 

In this section, the 2D TI method will be briefly reviewed first. The generation of the 3D 
mesh for the QQ/LC region with inclined slots is then given. 
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The Transfinite Interpolation Method 

Following Gordon and Hall [8], consider a 2D region as shown in Figure 3 represented by a 
vector- valued function F which is covered by a normalized parametric coordinate system in % and 
tj. Accordingly, the corresponding transformation plane (^-T| plane) is a square region. There are 
two sets of constraint curves embedded in ~F. One set has M constraint curves (solid lines) denoted 
by F&, Ti), where k = 1, 2,...., M, ^ = 0, and = 1. The other set has N constraint curves (dashed 
lines) and is denoted by F(£, T\j), where j = 1, 2,...., N, tj, = 0, and "n N = 1. The Boolean sum 
projector or transfinite interpolant, which matches F exactly along each of the M+N constraint 
curves, is constructed as follows [8]: 

{%®¥j[F] = %m + ¥jF]-t ?&[¥] (!) 

where the transformation projectors, and are defined in terms of blending functions 4> and y 
as 


^[Fj = S UWM, 

0 <rj < 1 

(2) 

^[F] = £ Vj(T\jF(£,T\j), 

;' = i 

0<£<1 

(3) 


and the product projection which matches F at the intersections of the two sets of constraint curves 
is defined by 

2VP„[F] = i . (4) 

In the above equations, the blending functions, <j> and \|r, are expressed as: 
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k = l,2,...,Af 


(5) 


♦*<&=4r 

rm-oj 

m = 1 
m*k 


ft 01 -'*'») 

Y/iD-TT* 

ncio-Y.) 

n = 1 
n*j 




where the coefficients <t> and T* are defined by 


(6) 


•k-£i=^T (7) 

m = \kmax - 1 

\p _ ^ ./(^i") j; (8) 

’ n = i jmax - 1 

Here kmax and jmax are the total number of grid lines in the respective ^ and T| directions. k(^ m ) 
and j(T| n ) are functions giving the k and j indices of the constraint curves, respectively. 

Equation (1) approximates F in the sense that it describes an approximate domain matching 

F on the non-denumerable (transfinite) set of points comprising the M+N constraint curves [8]. 
A pplication to the Combustors wi th Inclined Slots 

To facilitate the grid generation scheme using the TI method described above, the QQ/LC 
computational region shown in Figure 2 is divided into three subdomains: the RB zone together 
with the convergence section, the QQ zone with slot, and the divergence section along with the LC 
zone. The TI method is first applied to the QQ zone with grid lines exactly describing the desired 
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slot geometry, but with zero inclination angle. Next, the TI method is applied to the other two 
subdomains but with grids generated only for one r-z plane. A patching procedure is then used to 
generate the grid lines in the azimuthal (0) direction, based on the previously generated QQ zone 
3D grid. The final step is to twist the mesh by the desired angle, i.e., the slanted slot angle, to form 
the inclined slot configuration. Following is a more detailed description of the grid generation 
procedure. 

Referring to Figure 4(a), first, a rectangular "z-p" plane is formed by linearly connecting two 
points at the beginning and the end of the QQ zone. The slot is then projected from the curved 
combustor wall onto this z-p plane. Figure 4(b) shows the outcome of these steps. The parametric 
equations describing the four boundaries of the z-p plane shown in Figure 4(b) are 


z (£ rij) = (z f - z t )% +z if p (i;. Hi) = Pi 

z(^ , n 4 )=(z / -z i )^+z,, p(^)=p f (10) 

z(^,Ti) = z 1 ., p(Z >l ,r\) = (Pf-p i yn + Pi ( n ) 

z(&,r\) = z f , p&,r\) = (Pf-Pi)T\+P>' ( 12 ) 

The equations for the lines describing the slot boundaries are 

z(£,ti 2 ) = (z / -z,)^ + z J - , P&t\ 2 ) = Pc- q - 5S w (13) 

z&ti 3 ) = (z / -z i )^ + z ; , p(£ > ,T\ i ) = p c + 0.5S w (14) 

z(^ 2 ,Tl) = Z c - 0 . 55 i , P(^ 2 ,T\) = (Pf-Pi)T\+Pi ( 15 ) 

z(^3,ti ) = z £ + 0.5S l , p(^,T\) = iP/-Piyn + Pi ( 16 ) 


where z c and p c are the slot center coordinates and S w and S L are the modified slot width (WO and 
length (L), respectively, given by 
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(17) 


S w — 


w 

cos e slol 


S L — L ■ cos 6 s io t 


Here 0 slo , is the slot inclination angle. £ 
coordinates, 0 < T) < 1, and defined by 


$- 


T\ = 


Z/-Z; ’ 

P-Pi 

Pf-Pi 


z,<z< z, 


Pi <p <p f 


and Tl in Eqns (9)-(16) are the normalized parametric 


(18) 


These are the two sets of the constraint curves (M = 4 and N = 4) used in Eq (1) to generate the 
mesh shown in Figure 5, i.e., Eqns (11), (12), (15), and (16) form one set of the M constraint curves 
and, Eqns (9), (10), (13), and (14) form the other set. The input data for generating Figure 5 are 
given in Table 1 in which G seclor is the sector angle. Note that a cubic polynomial stretching function 
[13] is used to space the grid points in the \ direction so that the grid spacing gradually changes 

away from the slot as shown in Figure 5. 

Table 1 : Grid Data for Figure 5 


kmax = 22 

WW = ‘ 

II 

u> 

k&)=16 

k&) = 22 

o 

II 

N 

z r = D 

jmax = 19 

i(T|i) = 1 

i(Tlz) = 8 

i(n 3 ) = 12 

i(n 4 ) = 19 

Pi = 0 

Pr= [y(l - cos Qsector)] 1/2 
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With the z-p plane grid defined, an inverse procedure is carried out to project grid points back 
to the combustor wall (See Figure 6). Accordingly, the x and y coordinates of the grid points on 


the combustor wall must satisfy the following two equations simultaneously 


^ imax.j.k X imaz,j,k 


(D V 
V 2 ) 


yimax,j,k ^1 ^'’2p^imax,j,k 


1 < j < jmax ; 1 < k < kmax 


(19) 


where imax is the i index on the combustor wall and equals the total number of radial grid points. 
The coefficients C t and Q, in Eq (19) for each point are determined such that the projection line is 
normal to the z-p plane at the point of intersection. The z-coordinate on the combustor wall is 
simply the z-coordinate on the z-p plane. From Eq (19), 0 for each point on the combustor wall is 


determined by: 


0 . = tan 

v wrtax,y t * 


-1 


yimax ,j,k 
k J 


( 20 ) 


To obtain the x-, y-, and z-coordinates of grid points at the other radial planes, i.e., 0 < r < D/2, the 
following equations are used 


i,j,k r‘i,j,k^QS^vnaxJ,k 
yi,j,k = r i,j,k S in GinaiJ.* 


^ i,j,k ^ imax,j,k 


( 21 ) 


where ry k is calculated as follows: 
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1 < i < imax 


( 22 ) 


D_ 

^ imax,j,k ^ 

nj,k = n + ij,t- Ar i 

Here is the grid spacing in the radial direction. 

This completes the description of the grid generation for the QQ zone with zero slot inclination 
angle. In the case of non-rectangular slot configurations, e.g., circular holes, Eqns (13)-(16) need 
to be modified accordingly. 

To generate grids for the RB/convergence section, one r-z plane with 6 equal to 0 is used. 

Figure 7(a) shows the location and geometry of the appropriate embedded constraint curves, M of 
3 and N of 2 in this case. The TI method is then used to generate the mesh for this plane with one 
grid line exactly following the combustor wall in the convergence section. To obtain the grid for 
0 not equal to 0, a patching procedure is used to match the previously generated QQ zone 3D grid. 
Similar procedures are carried out for the divergence/LC section. Figure 7 (b) shows the embedded 
constraint curves (M = 3, N = 2) for this section. 

A perspective view of the grid at this point in the analysis is shown in Figure 8. Figure 9 
shows typical r-0 andr-z planar views of the mesh in Figure 8. There were 2 1 points in the r-direction, 
19 points in the 0-direction, and 53 points in the z-direction. 

To modify the grid to account for a slanted slot, the mesh is twisted by the desired slot 
inclination angle, 0 s iot, about the centerline and a reference base plane (r-0 plane). The following 
transformation equations then result: 
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XneW i.j.k = X i,j, 

k cosQ k -y it j 

,k sin©* 



y new ij,k = x ij 

,* sin ©*+?;,;, 

,k cos ©* 



ZneW i,j,k ~ Z iJ, 

k 


1 < i < imax ; \<j <jmax 

(23) 


where © k is the k* plane twisting angle, relative to the (k-l) th plane, and is defined by: 




^ imax,jmax,k ^ imax,jmax,k - 1 

'x x iirua,jmax,k "b yimaxjmax.k 


kbase <k < kmax 


Here kbase is the k index of the base plane which is a fixed plane, i.e., © kb ase - 0- 

Figure 10 shows aperspective view of the grid generated with 0 slo , of 45° in which the reference 
base plane is taken as the inlet plane of the QQ zone (i.e. the last plane of the convergence section). 


APPLICATIONS 

To see the effect of the twisted mesh on the solutions, a pure pipe flow with uniform inlet 
and without swirl was analyzed for the gnd system with and without grid twisting. It was found 
that the nonphysical swirl induced by the twisted mesh was less than 1% and, thus, could be 
considered insignificant. The approach was therefore considered appropriate for application to 
more realistic problems. The mesh shown in Figure 10 was used for the numerical solutions to be 
presented next. 

A 3D analysis of the reacting flow and heat transfer inside the QQ/ LC zones was studied 
using a modified version of the KIVA-II code and the grid-generation scheme developed here. Only 
a brief description is given here to facilitate the interpretation of the results; however, complete 
details can be found in [5-7]. 
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In the problem considered, a hot rich mixture from the RB zone enters the QQ section with 
inlet conditions obtained from the study described in [5,6]. There are 8 equally-spaced 45 slanted 
slots with slot aspect ratio (L/W) of 12 around the QQ zone, resulting in 9^ of 45°. The 
jet-to-mainstream mass flow rate and momentum flux ratios are 3 and 60, respectively, and the slot 
discharge coefficient is 0.6. At the slot, the dilution jet is air (consisting of 76.8% of N 2 and 23.2% 
of Oj) which enters through the slanted slot with uniform radial velocity. The following jet inflow 

conditions are used: 

Temperature = 1000 °F (811 °K) 

Pressure = 90 psia (6.205x1 0 6 dynes/cm 2 ) 

Turbulence intensity (v’) = 0.1 of V 
Turbulence length scale = 0.13 of D 
where V is the jet inflow radial mean speed. 

Figure 1 1 shows the computed isotherms and velocity vectors at the inlet plane which were 
taken from the study of the RB zone calculation [5,6]. Isotherms and velocity vectors at the center 
plane of the slot are given in Figure 12. Due to the slanted slot, neither the two vortices which form 
on both sides of the jet nor the isotherms are symmetrical with respect to the jet. It can be seen that 
near the wall the high velocity dilution jet deflects the main flow, causing swirl rotating in the 
direction parallel to the slot. However, near the axis, the main flow deflects the dilution jet, setting 
up a counter-rotating swirl in the direction of the main flow. 

The interaction of the inclined jet with the m ain flow is the cause of the co- and counter- rota ti n g 
bulk flow patterns shown in Figure 13 (isotherms are also given). This is the plane located about 
0.991 inches (2.517 cm) downstream of the slot center. From Figure 13, we see there is a high 
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shear layer (HSL) located near the mid-radius of the sector. The isotherms plot indicates that the 
HSL mainly contains lower temperature fluid, an indication of the effect of the jet penetration. The 
sandwiched-ring-shape (for a whole circle) temperature field, in the order of hot-cold-hot, can also 
be seen clearly in Figure 1 3. This bulk swirl flow and the sandwiched-ring-shape isotherm patterns, 
which are typical of confined slanted jet-in-cross flow, have been observed experimentally [14]. 

The bulk swirl flow persists further downstream, as shown in Figures 14 and 15. From these 
two figures, one can see that the HSL is moving outward radially before entering the divergence 
section, see Figure 14. However, at the end of the divergence section (Figure 15), the HSL is shifting 
toward the axis (radially inward). 

At the exit of the LC zone, the bulk flow has almost disappeared and the fluid temperature 
is in the range of 1060 to 1340 K, see Figure 16. A simple one-dimensional analysis (Figure 17) 
will show that, with the given input conditions, the outlet equilibrium temperature should be 
approximately 1208 K. 

CONCLUSIONS 

We have presented an efficient algebraic grid generation scheme based on a 2D transfinite 
interpolation technique to generate the 3D grid for gas turbine combustors with slanted slots. This 
scheme can easily be extended to gas turbine combustors with slot or other dilution hole geometries. 

To demonstrate the feasibility and the usefulness of the technique, it was applied to the analysis 
of the QQ/LC zones of an STC combustor. It was found that the fictitious swirl induced by the 
twisted mesh was less than 1% and, thus, could be considered insignificant. Major features of the 
flow and temperature fields taking place inside the QQ/LC zones were satisfactorily modeled. 
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Formation of the co- and counter-rotating bulk flow patterns and the sandwiched-ring-shape 
temperature field that occurs in the model can be seen clearly, which is consistent with experimental 
observations. 

Numerical solutions show that the method presented can be used as a tool to generate the 
grid points for gas turbine combustors with slanted slots. 
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QQ/LC Zones 

FWRB Zone i 



Figure 1 . — Schematic of a staged turbine Combustors. 



Section A A 

Figure 2.— Geometry of the system under investigation and its 
specifications. 



Figure 3.— Diagram of the embedded constraint curves in the parametric cordinate system. 



Figure 4.— The z-p projection plane and the embedded constraint curves showing the location 
and size of the slot. 
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Figure 7. — Diagrams of the RB/convergence sections at the Figure 8. — Prospect view of the grid system before twisting. 

divergence/LC section constraint curves. 
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Figure 1 0. — Perspective view of the grid system after twisting. 
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Figure 15.— Velocity vectors and isotherms at the end of the diver- 
gence section. 
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Figure 1 6.— Velocity vectors and isotherms at the outlet 
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Figure 17.— A control volume used to estimate the equilibrium temperature by energy 
balance. 
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